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We present 90% confidence level limits on magnetic monopole production at the Fermilab Tevatron 
from three sets of samples obtained from the DO and CDF detectors each exposed to a proton- 
antiproton luminosity of ^ 175 pb~^ (experiment E-882). Limits are obtained for the production 
cross-sections and masses for low-mass accelerator-produced pointlike Dirac monopoles trapped and 
bound in material surrounding the DO and CDF collision regions. In the absence of a complete 
quantum field theory of magnetic charge, we estimate these limits on the basis of a Drell-Yan 
model. These results (for magnetic charge values of 1, 2, 3, and 6 times the minimum Dirac charge) 
extend and improve previously published bounds. 
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I. INTRODUCTION 

The existence of even one magnetic monopole with 
magnetic charge g explains the quantization of electric 
charge e in terms of the Dirac quantization condition 1] 
eg = nhc/2^ n = ±1, ±2, .... (Throughout this paper we 
use Gaussian units, but numerical results are expressed 
in SI.) Besides explaining the quantization of electric 
charge, the existence of magnetic charge results in the 
dual symmetrization of Maxwell's equations and is 
not forbidden by any known principles of physics. The 
minimum magnitude of the quantization number is n = 1 
according to Dirac orn = 2 according to Schwinger 3] . If 
e is the charge of the electron, these magnitudes become 
n = 3,6, respectively, if quantization via quark electric 
charges is possible. It should be emphasized that mag- 
netic charge, like electric charge, is absolutely conserved, 
so the lightest magnetically charged particle is stable, 
unless annihilated by its antiparticle. 

Throughout this paper we refer to magnetically 
charged particles as magnetic monopoles, or simply 
monopoles. However, as Schwinger emphasized 0, 
magnetically charged particles could also carry electric 
charge; such particles he christened dyons. The quanti- 
zation condition for a pair of dyons labeled 1 and 2 is 
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(1.1) 



We will not explicitly mention dyons further. We merely 
note that the considerations here should supply similar 
cross section and mass limits on dyons as for monopoles. 
There would be modifications for dyons — for example. 
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there would be a binding contribution due to electric 
Coulombic attraction or repulsion, but because electric 
charges are so much smaller that magnetic ones, the mag- 
netic contributions are overwhelming. More significant 
probably are acceptance modifications, due to changes 
in energy loss, but we expect the quantitative impact of 
these changes to be small (less than 10%.) 

If they exist, monopoles will presumably be abundant, 
or can be pair-produced by some appropriate mechanism, 
and be trapped in matter. Previous (direct) searches for 
trapped and bound magnetic monopoles in various accel- 
erator samples H H 0, Q, H, |1M El E H III El , in 
meteorites 16j, and lunar soil UJ, as well as an earlier re- 
sult [l^[l^ from this experiment have been made. Other 
(indirect) searches by other methods are not covered in 
this paper, but are reviewed elsewhere II^ [201 . Here we 
report 1) a reanalysis of the data of Ref. [13, 2) the data 
of Luo [13 as well as 3) that of a third set of samples from 
the CDF detector recently measured and analyzed. This 
extension of limits is experimentally driven. Theoretical 
motivations derive from the expectation that monopoles 
from spontaneous electroweak scale symmetry breaking 
might give rise to monopoles of mass ~ 2.5 to ~ 15 TeV 
21, 22], although we here can only search out to a mass, 
in our Drell-Yan modeling, ~ 0.4 TeV. The LHC using 
this method would allow one to approach 2 TeV. Unfor- 
tunately no accelerator is currently envisioned that will 
reach the theoretically interesting region of 10 to 15 TeV. 
In view of our nearly complete absence of knowledge of 
the origin of particle masses, we should not exclude any 
mass region from an experimental search. 

This paper presents the experiment, analyses, and re- 
sults of our search for monopoles. Section ^ outlines 
the basics of the detector apparatus; details are provided 
elsewhere 19j. Section UTTI covers the calibration and lin- 
earity of the detector. Section HVl discusses the samples 
and monopole energy loss leading to the stopping and 
capture of monopoles in the sample material. Section 
IVl describes the analysis of the measurements. Section 
[^discusses the transformation of the data to monopole 
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FIG. 1: Schematic radial cross-section of the monopole de- 
tector at a 2:1 width to height ratio. The elements generally 
are rings, tubes, or cylinders concentric with the indicated 
center line. The height, which was 155 cm during the set 1 
measurements of the DO aluminum samples, was increased to 
185 cm for the subsequent CDF lead and aluminum samples. 



cross-sections and mass limits. Finally we summarize in 
Section EUl 



II. MAGNETIC MONOPOLE DETECTOR 



We use the induction method of Alvarez, et al. j23ll24| 
to detect monopoles. A large warm bore cryogenic de- 
tector, similar to that of Jeon and Longo 16], was con- 
structed at the University of Oklahoma. The active el- 
ements of the detector, shown schematically in Fig. ^ 
are two 19 cm diameter superconducting loops each con- 
nected to a dc SQUID (superconducting quantum inter- 
ference device). The magnetic flux from a magnetic mul- 
tipole passing through the loop induces a change in the 
loop's supercurrent because the Meissner effect prevents 
a change in the magnetic flux through such a loop. If the 
multipole is dipole or higher, then the total net change 
when the sample traverses the loop starting and ending 
at large distances goes to zero, although locally super- 



currents are present when the sample is in the neigh- 
borhood of the loop. But if a monopole passes through 
the loop, then a nonzero net change in the supercurrent 
occurs, giving a current "step" characteristic of a mag- 
netic charge. This effect can be observed from even one 
monopole in a macroscopic sample due to the long range 
(~ 1/r^) nature of the associated magnetic field. The 
change in supercurrent is detected by the SQUID con- 
nected to the loop and converted to a voltage by the 
squid's preamplifier and controller (see below). 

Measurement of samples of a size less than 7.5 cm in di- 
ameter by 8.5 cm in length is made by repeatedly passing 
them through the 10 cm diameter warm bore centered on 
and perpendicular to the loops. A vertical excursion of 
1.1m around the position of the loops is typically made. 
In a central 65 cm region this allows for the magnetic 
effects of induced and permanent dipole moments of the 
sample to start and return to zero on each up and each 
down traversal taking some 25 s each. A net data rate of 
10 Hz is recorded for each of the SQUIDs. Also recorded 
are the readings of an accelerometer, the vertical position 
of the sample recorded by an optical encoder, the num- 
ber of increments taken by the stepper motor moving 
the sample, and the time. The data acquisition (DAQ) 
was performed with Apple Macintosh computers running 
under National Instruments LAB VIEW programs 25]. 

Other magnetic and electronic signals affect the 
SQUIDs also, causing systematic errors. One has to deal 
with: 

1. Permanent dipoles in the samples (presumably mi- 
croscopic particles of magnetite or other ferrites) 

2. Induced dipoles (because the samples are conduct- 
ing metals) 

3. Contamination of the transporting nylon string and 
copper wire 

4. Ground loops 

5. External electronic device interference (certain 
clocks, welding operations, etc.) 

6. Thunderstorms (which forced suspension of opera- 
tions) 

7. Mechanical vibrations (external limited by 
dampers and isolation, internal due to cryogen 
boiling) 

8. Small variations in the warm bore magnetic field 
gradient 

9. Unidentified sources (some days operations had to 
be suspended). 

In addition to data running, background running was 
required, which was subtracted from the data runs. Even 
during good running conditions, source (3) was always 
present to some degree. Source (1) was present for most 
samples; large dipole moments off-center in the volume of 
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the sample give dipole tails that are unbalanced and can 
mimic monopole steps. Samples with large dipole sig- 
nals need to be vetoed, because the SQUID loses count 
of the number of flux quanta. Induced signals, source 
(2), which are oppositely directed on up/down traver- 
sals, are pairwise cancelled. The temporal dependence of 
these signals [a complete sample measurement takes some 
twenty (20) minutes] is minimized by using only time- 
coincident traversals from the two squids taking data. 
It is to be noted that source (1) is variable. On re- 
peat traversals, the magnitude of the dipole signal may 
change, most likely because the magnetite grains may be 
re- magnetized. The induced dipole signal, source (2), is 
proportional to the conductivity of the sample, the mag- 
netic field gradient in the warm bore, and the velocity of 
the sample. The velocity used (one meter per 25 seconds) 
kept the magnitude smaller than the permanent dipole 
value in general. 

Sporadic electrical and electronic interferences listed 
above occasionally caused automatic resets of the SQUID 
controllers. The operator had to restore the output levels 
of the SQUIDs manually and continue the measurement. 
Later, in the analysis stage, these reset (bad) data traces 
are vetoed pairwise. 



sensor, (i^squid is the fiux change in the SQUID, (i^ioop is 
the flux change in the detection loop, 0o is hc/2e = 2.07 x 
10~^^ Wb, the superconducting fluxoid unit, and dV is 
the change in output voltage of the SQUID controller. 
Our SQUIDs were two junction DC types manufactured 
by Quantum Design [2^ . 

dV/d^squid is the so-called transfer function, "XF," 
measured by a test pulse injection procedure given by 
the manufacturer, for our SQUIDs approximately 0.5-0.8 
V/^o- <^^squid/<^^ is a parameter measured and furnished 
by the manufacturer, dl/d^ioop is calculated as follows. 
Since the magnetic flux linking a superconducting loop is 
conserved as a result of the Meissner effect, we have by 
integrating the generalized Maxwell equation 
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(3.2) 



Jm being the magnetic current density, the following re- 
lation for the current induced in the detection loop by 
one flux quantum passing through it (see Appendix): 
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III. CALIBRATION AND LINEARITY 

The schematic circuitry of the SQUIDs and the asso- 
ciated detection loops are shown in Fig. [21 During op- 
eration of the experiment, the SQUIDs are tuned and 
their transfer functions measured periodically according 
to the manufacturer's speciflcations to keep them oper- 
ating with constant sensitivity. The sensitivity can be 
approximately predicted from the SQUID sensor's pa- 
rameters and accurately calibrated by measurement. The 
absolute calibration of an expected signal from a Dirac 
monopole is made using a "pseudopole." A long thin 
magnetic solenoid carrying a small known current gives 
a calculable pseudopole (ps) at either end. The pseu- 
dopole can either be passed through the warm bore of 
the detector in a way similar to the samples (discussed 
later), or it can be placed in a given position with one 
end fully extended through the SQUID loops and the 
solenoid current repeatedly switched on and off. Both 
methods were used. We now discuss the parameters and 
the calibration. 



The expected response function (RF 
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the SQUID system depends on the following parameters 
(where the current to voltage conversion is done by the 
SQUID controller's ampliflers): 

dV d^squid dl . 
^^ = 3^^ 77 — 1^ V^o, (3.1) 

where dl is the change in current induced in the SQUID 



where (see Fig. |2| Li is the inductance of the detection 
loop, L2 is the inductance of connecting twisted pair, L3 
is the input inductance of the squid sensor, r is the radius 
of detection loop, and a is the radius of the superconduct- 
ing lead (Pb) shield. The factor (1 — {r/a)'^) inserted in 
Eq. (|3.3|) corrects for the trapping of flux coupled into 
the detection loop due to the cylindrical superconduct- 
ing shield. (See Appendix.) 

The sensor's input coil inductance L3 is measured to be 
of order 1.8 /iH (see Tabled). The inductance L2 is zero, 
and Li (of order 1 /iH) is calculated from the standard 
formula Mm : 
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where r is the radius of the loop and p is the wire radius. 
For a uniform current density in the loop C = 7/4; for a 
(superconducting) surface current C = 2 (a reduction in 
the inductance of 4 percent). In addition, the supercon- 
ducting shield reduces the value of Lioop by 9.5 percent: 



Li = Lioop(2)(l - 0.095) = 0.70/iH. 



(3.5) 



This comes about as follows: The inductance between the 
detection loop and the superconducting shield is given 
by the following formula, rather easily derived from the 
inductance between two loops [28^ : 
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FIG. 2: Schematic of the monopole detector's signal measuring electronics. The flux passed through the superconducting loop, 
LI, is converted into a current change that induces a flux change in the DC SQUID ring. An induced feedback current is 
amplified and converted into the output signal voltage in the preamplifier/controller and processed by the LABVIEW data 
acquisition (DAQ) program running in the Macintosh computer. 
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where a is the radius of the shield and 21 its length. In- 
verting the inductance matrix changes the effective in- 
ductance of the detection loop by a significant amount 
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and putting in the values a = 14.73 cm, r = 9.855 cm, 
/ = 30 cm, and p = 0.02 cm gives Li as in Eq. (|3.5|) . 
where Lghieid = 47ra^/(2/c^) is the self-inductance of the 
superconducting cylindrical shield. 

The experiment had four (4) SQUID sensors referred 
to as DCl, DC2, DC3, and DC4. DCl failed at some 
early time, was repaired and subsequently called DCIR. 
Most of the sample measurements for the final analysis 
were performed with the pair DCIR and DC2. However, 
some measurements were made with other pairings of the 
four SQUIDs, and in a number of configurations, which 
are catalogued in the tables which follow. In initial tests, 
we had an emulation setup with a small bore and small 
detection loops immersed in liquid helium contained in a 
"Research Dewar." We also had two current sources for 
the pseudopole, one having multiples of a current unit [~ 
1.55 V battery/500 M^l, equal to about 0.7 Dirac poles — 
see Eq. (|3.9|) below] equal to 1, 2, 5, 25, 100 current units 
and a later one with multiples of 1.7, 2.5, 5, 50 current 
units. We also had an early pseudopole of 0.5 meter 
length and a later one of 1.016 meter length. The actual 
currents for each number of units was measured with a 



picoammeter to be correlated with the voltage response 
of the calibration measurements. 

The individual measured transfer functions, XF, in 
V/(j)o are given in Fig. El We see two lines of values at 
0.75 and 0.64 for most of the plot, with values of 1.5 and 
1.3 at year 2002.5, where we doubled the output gain of 
the squid controllers. The double valuedness is due to the 
presence or absence of an external passive RC low pass 
filter, which exhibits attenuation to the signal before be- 
ing converted on the ADC board in the computer. (The 
controller has an internal active low pass filter which can 
be switched on or off, which in comparison is attenuation- 
free well below the rolloff point.) Tabled identifies in the 
name DCi the filter configuration as "NF" or "EE" for 
"no filter" (or non-attenuating active internal filter) or 
"external filter," respectively. In the table, *2 or RD in- 
dicate the doubled gain or early Research Dewar data. 
A gain 0.5 data point was taken while testing modifica- 
tions to the controllers, and is exactly half of the gain 
1.0 point above it (see Fig. The transfer functions 
and their errors shown in Fig. |31 have been averaged as 
appropriate and entered into Tabled Tabled also shows 
the squid /dl values, the input coil inductances L3, the 
dl /d^ioop values calculated from Eq. (|3.3|) . and the theo- 
retical RF's [Eq. (|3.1|) ] shown (mV/Dirac pole). Note one 
Dirac pole corresponds to two fluxoid units, Ang = hc/e. 

The experimental data (mV/nA) for traversals of the 
pseudopole through the warm bore, as with samples, and 
for the pseudopole "parked" at one position and its cur- 
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FIG. 3: The myriad of transfer function (XF = dV/d^sciuid) measurements made over the course of the monopole experiment. 
Time intervals for the SQUID configurations in use at the various times, as well as the periods of measurement of the different 
sample sets are also indicated. RD = research dewar, DO and CDF are the collider experiments from which the samples were 
taken, RM = remeasurements and 1, 2, 3, 4 refer to the various SQUID probes (see text). 



TABLE I: Transfer functions and predictions. The subscript s denotes SQUID. Here {Li, 1 - (r/a)^} = {0.70 /iH, 0.552} for 
all entries except for the first (RD) where they are {0.01 /iH,0.87}. Remember (I)d = 2$o- A{dV/d^s) is the error in dV/d^s, 
and likewise A (RFth) is the error in RFth- 
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0.12 



rent turned on and off, are shown in Fig. ^ and 
respectively. The two methods agree within experimen- 
tal uncertainties. One again sees (at early times) the 
RD data, then the normal mid- year data, and finally the 
doubled gain data in 2002. Averaging the appropriate 



groupings of DCi gives the values (mV/nA) in Table ITTI 

The conversion from the response in mV/nA to that 
in mV/Dirac pole is given by considering the pseudopole 
solenoid's magnetic moment both as a current loop and 
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TABLE II: Measured pseudopole responses and comparisons to predictions (Table 0. (Note: mV/nA x nA/(j)D — mSf /(j)D-) 
Here "Ratio" is the ratio of the theoretical response RFth given in Table |l] to the measured response/pole. 
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for the new (old) pseudopole, respectively. 
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FIG. 4: Individual measured pseudopole responses (mV/nA) 
for a) traversing the loop similarly as the samples and b) 
turning the current on/off while the pseudopole is "parked" 
at one position. 



as a pole-antipole dipole, i.e.. 



1 



gL = -NiA 

c 



(3.8) 



where the length of solenoid is L = 1.016 (0.50) m for the 
new (old) pseudopole, respectively, Ni is the ampere- 
turns of current, N = 4710 (2440) turns, and the cross- 
sectional area of the solenoid is A = 1.533 ± 0.011 x 



Using these conversion values we get the measured 
pseudopole response in mV/^D pseudopole values ver- 
sus DCi configuration shown in Table ITTI and Fig.ISK and 
the expected/measured ratio in Table Hll and Fig.ISt), re- 
spectively. A weighted average of these ratios for the 
nine values is 1.21 ±0.07. The discrepancy is presumably 
due to approximate values of Quantum Design's mea- 
sured parameters, and the elementary treatment exclud- 
ing shielding internal to the SQUID probes, etc. The 
values plotted in Fig. [SJi are used to determine cuts on 
the possible monopole "steps" seen in the analysis of the 
data given later below. Data taken at year 1997.7 give 
the linearity response shown in Fig.[SJi andlHt). One sees 
that the calibration is very linear, and that signals down 
to one-half a Dirac pole can, in practice, be measured. 

The shape of an expected monopole step is obtained by 
a subtraction of two runs with different pseudopole cur- 
rents, as shown in Fig.[2K- Five such shapes are compared 
in Fig. 03. The theoretical modeling, given in the Ap- 
pendix, is compared to the experimental data in Fig.Q). 
The agreement between the theoretical and experimental 
shapes is good. We conclude that the calibration, lin- 
earity, and positional response (shape) indicate that the 
experiment is understood, and that operating conditions 
over six years are reasonably consistent and stable. 
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IV. SAMPLES, STOPPING AND TRAPPING 
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A. Samples 

There were three sets of samples obtained from dis- 
carded material from the upgrading of the DO 30] and 
CDF ^1] detectors: 1) Be beam pipe and Al "extension" 
cylinders from DO, 2) Pb from the forward/backward 
"FEM" (Forward ElectroMagnetic calorimeters) of CDF, 
and 3) half of the Al cylinder ("CTC" support) from 
CDF. Sample set 1 was initially reported 18], 2 was given 
in Luo's PhD. thesis 

'.19], 

and set 3 here. In this paper, 
all three sets are being reported according to a final con- 
sistent analysis. 

Sample set 1 comes from the two Al extension cylin- 
ders (extending beyond the main central detectors of DO 
inside the liquid argon calorimetry), each of 150 cm di- 
ameter by 46 cm length by 1.26 cm thickness. They came 
as 16 plates. Four of these plates were further cut by wa- 
ter jet, whereas the balance were cut on a bandsaw; no 
magnetic degradation was observed from the bandsaw- 
ing in comparison to the water jet samples. These cut 
pieces were 7 cm by 7.6 cm or 6 cm by 7.6 cm in size. 
Two of each were bundled into a "cylinder" constituting 
a "sample." There are a total of 222 Al samples. The Al 
sawings and other small scrap pieces are accounted for 
later. 

In addition, sample set 1 included the DO 5 cm diam- 
eter, 0.05 cm thickness Be beam pipe. The central 46 
cm section, centered on the collision region and covering 
nearly the full solid angle, was cut into six 7.6 cm pieces. 



Sample set 2 comes from twelve 2 m by 2 m by 0.5 cm 
thick Pb sections cut from the center of the full 3.04 m by 
3.04 m Pb layers. These twelve layers were those closest 
to the interaction region of which six were located on the 
east side and the other six on the west. Each layer had 
an octagonal hole at the center (of approximately 30 cm 
"diameter") for the passage of the p-p beams, etc. Each 
2 m by 2 m section came as six pieces of 0.667 m by 1 
m, labeled T (top), M (middle) or B (bottom) and also 
labeled N (north) or S (south). The six sections were 
also labeled as E (east) or W (west) and by layer number 
(1,2. .. ,6). These sections were sheared into twelve strips 
of 8.3 cm width across the 1.00 m dimension, and each 
strip was rolled into a cylinder of approximately 7.5 cm 
diameter. These constitute some 816 samples. A typical 
layer is shown in Fig. [HI 

Of the 816 samples, only 664 have been successfully 
measured. The rest were unmeasurable due to huge 
dipole moments caused by the (ferromagnetic) red paint 
on the first layer samples. Fortunately, layer 1 samples 
turn out not to be needed. The large solenoidal mag- 
netic field of CDF accelerates the monopoles to a high 
enough energy (70n GeV additional) that they penetrate 
the early layers, reaching layers 4, 5, 6 for n = 1 and layer 
2 for n = 2, 3, and not at all for n = 6. The scrap pieces 
and other unmeasurable samples of Pb are accounted for 
later. 

Sample set 3 comes from the CDF's inner support 
cylinder, which has eight 0.63 cm thick cover plates on 
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FIG. 7: a) Pseudopole response, SQUID DC2, versus position 
for two currents, ±50 current units (about ±35 Dirac poles) 
and their difference, giving the experimental monopole step 
shape, and b) several such shapes [one of DCIR, three of 
DC2 including 7a) above, and one of DCS] compared, and c) 
Fig. [7Ji compared to the theoretically expected shapes: loop 
in free space (poor) and in a superconducting can (SC) 
(good)-see Appendix. 



a cylinder of 2.74 m diameter by 2.64 m length. Six 
of these cover plates, covering three-quarters of the az- 
imuthal angle about the beam direction, were made into 
samples, similar to the DO ones. These were sheared on 
a large shearing machine, first into long strips and then 
each strip into shorter sections. Eight pieces (four of each 
size) were bundled into a 7.5 cm diameter by 7.6 cm long 
cylinder (the CDF Al being half the thickness of the DO 
extension pieces) yielding 404 samples with little scrap. 
Due to limitations on time, funding and personnel, 132 
of the samples did not get measured, leading to only one- 
half the azimuthal solid angle being covered. These and 
the scrap are accounted for later below. 

It turns out that all of the CDF Al samples were very 



M 



B 




FIG. 8: A lead sample layer is cut into 68 samples. The 
Monte Carlo simulation shows that only the samples (52 in 
one layer) subtending an angle 3° < < 7.5° have signifi- 
cant acceptance and should be considered. The cross-hatched 
samples and others with cuts on large dS (statistical) and sS 
(systematic) errors are excluded from the analysis; the (solid 
angle) coverage falls to 53% for n = 1 monopoles (see text). 

magnetic (similarly to the CDF Pb red painted layer 1 
samples). We had to demagnetize them. Degaussing 
with an AC field coil was ineffective. But as ferrites have 
Curie points below 585 °C and the melting point of 
Al is 660 °C (as opposed to 327 °C for Pb), we were 
able to demagnetize all but a few of the samples by a 
heat treatment, "soaking" the samples at 610 °C for one 
hour.^ The dipole amplitudes decreased from an unman- 
ageable 3-10 V to the usual 30-100 mV. 



B. Stopping of the Monopoles 

1. Energy Loss 

The energy losses {dE/dx) of magnetic monopoles 
traversing material absorbers are caused by the inter- 
action of the moving monopole charge {g = ngo) with 
the electric field of the atomic electrons, i.e., the x E 
Lorentz force. The velocity dependence of this force can- 
cels the 1/ (velocity)^ dependence of the usual charged 
particle dE/dx. Either classicall y 1271 . quantum mechan- 
ically Is^, or field theoretically |35j . approximately one 
simply substitutes {gPY for {zeY in the usual charged 
particle dE/dx formula. Kazama, Yang, and Goldhaber 



The resulting small increase in thermal energy is completely neg- 
ligible compared to the binding energy of monopoles to matter. 
Monopoles bound with a keV of more of energy should be per- 
manently trapped in the material. See Ref. |33i|| . 



9 



have obtained the diflFerential scattering cross sec- 
tion for an electron moving in the magnetic field of a 
fixed magnetic pole. Ahlen has used this cross- 

section to obtain the following expression for monopole 
stopping power: 



dE 
dx 



2J1 



-NA In 



B(\n\) 



2 2 



(4.1) 



where A^e is the number density of electrons, / is the 
mean ionization energy, K(\n\) = 0.406 (0.346) is the 
Kazama, Yang and Goldhaber correction for magnetic 
charge n = 1 (n > 2) respectively, b is the usual density 
correction and B(\n\) = 0.248 (0.672, 1.022, 1.685) is 
the Bloch correction for n = 1 (n = 2,3,6), respectively 
|3^ . (Of course, one must divide by the density to get 
dx in g/cm^.) This formula is good only for velocities 
/3 = i;/c ^ 0.1. For velocities /3 ^ 0.01, we use Eq. (60) 
of Ref. [43 as an approximation for all materials: 



dE 
dx 



(45 Gev/cm)n^/3, 



(4.2) 



which is linear in in this region. The two dE / dx veloc- 
ity regions are joined by an empirically fitted polynomial 
in the region of = 0.01-0.1 in order to have a smooth 
function of (3. For the elemental and composite mate- 
rials found in the DO and CDF detectors, we show the 
resulting dE/dx curves we used in Fig.jUl (See Ref. 19].) 



2. Stopping of Possible Monopoles 

The trajectories of possible monopoles are generated in 
a Monte Carlo program from the collision point through 
the appropriate detector elements slowing to a stop in the 
sample layers, including the acceleration (or deceleration) 
of the monopole along the external magnetic field lines of 
the CDF detector. The average polar angle relative to 
the beam direction is different for the three sample sets: 
= 35°, 6.5°, and 90° for sets 1, 2, and 3, respectively. 
The amount of material along a "typical" trajectory is 
shown in Table Iml for DO and Table ITVl for CDF. 



3. Trapping of Monopoles 

The monopoles having been stopped can then bind 
to the magnetic moments of the nuclei of the material 
present, and are therefore trapped. The interaction of 
the monopoles with the magnetic moments of nuclei and 
electrons can be strong enough to produce bound states 
under certain conditions and be trapped, having a very 
long lifetime in such bound states j^^- The trapping 
efficiency is very high. The theoretical modeling that 
has been done assumes "rigid" extended nuclei with or 
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FIG. 9: dE/dx vs P v/c) for a) elements and b) compounds 
and mixtures involved in the energy losses of monopoles pass- 
ing through the various detector materials. To increase legi- 
bility of the elemental figure, we have not shown the energy 
loss for Si and Be, which are similar to Al, nor for N2, which 
is similar to C. For the compounds and mixtures, we note 
that C3H6, CO2-C2H6-AI, and Rohacell are similar to Hon- 
eycomb, and Si02 is similar to AI2O3. These dE/dx curves 
are for a magnetic charge value of n = 1; apart from the cor- 



rection terms K{\n\) and B{\n\ 
magnetic charge values. 



we multiply by n for larger 



without repulsive barriers, and some relativist ic calcu- 
lations have also been carried out. Electrons can bind 
to monopoles in a total energy zero state; this probably 
produces a small mobile system which will transfer to a 
nuclear magnetic moment leaving it bound to a fixed nu- 
clear site and permanently trapped. Thus we assume all 
monopoles bind to appropriate nuclei, i.e., those whose 
nuclear gyromagnetic ratio is sufficiently large (anoma- 
lous). These models predict, as summarized in Ref. |33( . 
that binding should occur for 2^A1 (100% natural abun- 
dance) and 2°^Pb (22% natural) but not for ^Be (100% 
natural). However, the estimated binding energies, e.g., 
0.5-2.5 MeV for aluminum, are large and comparable to 
shell model splittings, so we believe that in the presence 
of the monopole the nucleus will undergo nuclear rear- 
rangement and binding should in general result, even for 
^Be. Even an unreasonably small estimate for the bind- 
ing energy of 1 eV would give a lifetime of 10 yr j3^ . 
We have, therefore, good reason to believe that stopped 
monopoles will be trapped by the magnetic moments of 
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TABLE III: The materials and path lengths at an angle of 35° relative to the beam direction for the DO aluminum sample 
measurements (set 1). 



Detector Material Path lengths (g/cm^) 

Be pipe 0.16 

VTX C2H6 and CO2 gas, Al wire mixture 0.12 

TRD Al window 2.18 

N2 gas 0.01 

Polypropylene 1.12 

Xe, CH4, C2H6 gases 0.02 

Honeycomb 1.99 

CDC Ar, CO2, CH4 gases 0.07 

Kapton 0.1 

Rohacell 0.49 

Endplate (Al) 4.18 

GIO 1.11 



Total including Sample layer (Al) 11.6-17.6 



TABLE IV: The materials and path lengths in the CDF detector, at an angle of 6.5° from the beam for CDF lead samples 
(set 2), and at an angle of 90° from the beam for the CDF aluminum samples (set 3). The Pb samples are spaced by 0.6 g/cm^ 
of drift chambers. 



Detector 

svx 



VTPC 
FTC/CTC 



Total including 



Material 
Be pipe 
C 
Al 
Si 
AI2O3 
C 
Al 

Pb layer and drift chamber 
Sample layers 



Set 2 

Path lengths (g/cm^) at 6.5° 
0.82 



0.29 
3.57 
FTC: 0.25 
6.5 



Pb layer 2 
Pb layer 4 
Pb layer 5 
Pb layer 6 



11.4- 17.2 

24.5- 30.2 
31.1-36.8 

37.6- 43.3 



Set 3 

Path lengths (g/cm^) at 90° 
0.092 
0.209 
0.135 
0.210 

0.312 
CTC: 0.293 



Al: 1.25-2.97 



nuclei. 



V. ANALYSIS OF THE DATA 

The data analysis proceeded generally as follows, with 
modifications as required for the different sets of sam- 
ples. The time sequences of the SQUID's outputs were 
examined interactively, and bad sections, primarily sec- 
tions that contained a SQUID reset (see Sec.^, deleted 
pairwise (to keep induced signals canceling out pairwise) 
and the corresponding traces on the other SQUID also 
deleted (to keep them temporally consistent); typically 
80-90% of the traversals remained. A pedestal value, the 
SQUID output near the top end of each traversal (bot- 
tom for sample set 3), was subtracted from every voltage 
value along that up (or down) traversal. The values for 
each of some 90 (75 for sample set 3) small ranges of 
vertical positions were averaged, removing most of the 
random drift of the SQUIDs. The two SQUIDs' data 



were averaged, shifting one relative to the other by 10.1 
cm in position in order to superimpose their dipole re- 
sponses. The background samples were analyzed simi- 
larly and local groups of background runs were averaged. 
These background runs were subtracted from the sam- 
ples' spectra. A pair of horizontal lines was fit to two 
regions, one at the lower position and one at the upper. 
The difference in values of these two flat fits gives the 
step for that sample. Examples of processed data are 
shown, one for each sample set, in Fig.^J The steps for 
each sample set are histogramed and analyzed for con- 
sistency with the null hypothesis, absence of monopoles. 
The background subtraction ensures that the distribu- 
tion of steps centers on zero, since various background 
effects, such as the small effect of the magnetized string 
holding the sample in set 1, has been removed. These 
analyses were performed with a large set of macro proce- 
dures written for the program WaveMetrics' IGOR-Pro 
I^H running on Apple Macintosh computers. 

For sample set 1 (DO aluminum and beryllium) the 
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FIG. 10: Typical step plots: a) set 1, DO aluminum, sample 
"S115," b) set 2, CDF lead, sample "WSM28," and c) set 3, 
CDF aluminum, sample "A012." The steps are —0.4, —1.1, 
and 2.0 mV, respectively. 



background sequence was one background sample (al- 
ways the same sample) between every two samples mea- 
sured. For sample set 2 (CDF lead) it was one back- 
ground between every three samples. And for sample set 
3, the other samples were used as backgrounds for any 
given sample. The background averaging was done over 
time blocks of data. Each time the measurements were 
interrupted for unmanned time, a liquid helium fill, a 
thunderstorm, etc., some possible change in the environ- 
mental background was likely, so these were considered 
time blocks that required an independent background 
subtraction. 

Sample set 2's sample holder was a threaded metal 
rod with hexagonal nuts clamping the rolled lead sample 
cylinder. The overall magnetic dipole background was 
generally dominated by the sample holder, making the 
analysis of the CDF lead samples a bit harder than those 
of sample set 1. 

Sample set 3 differed most in its procedures. A greater 
signal sensitivity should give a better overall signal-to- 
noise ratio for each SQUID's data. As a major part of 
the noise was external to the SQUID electronics, the over- 




step (mV) 

FIG. 11: Histogram of the steps from 222 aluminum and 6 
beryllium samples from the DO experiment. The total of 228 
samples, with an rms spread of 0.73 mV, is compared to a 
Gaussian with that same standard deviation. 



all gain of the controllers was increased (see Section ITTTI 
above) by doubling feedback resistor values and halving 
capacitor values of the final operational amplifiers. This 
sensitivity gain apparently would have occurred. How- 
ever, as the experiment was winding down, with time, 
funding and manpower coming to an end, the failure of 
one of the two SQUIDs for set 3 was not repaired (in con- 
sideration of the time that would have been lost to warm- 
up, repair, cool-down, etc.). The remaining SQUID had 
a small 0.4 Hz oscillation which we then dealt with in 
the analysis software, by smoothing the spectrum over 
a number of the 10 Hz data input values before edit- 
ing the tr aver sals. In addition, we extended the range 
of the tr aver sals to 1.2 meter, at a larger stepper motor 
interval, leading to fewer (the 75 bins mentioned above) 
vertical position bins. The IGOR macros were modified 
to allow the computer to delete pairs of traces based on 
a chi-squared test of the consistency of differing traces 
in a sample's traversals from the average of those mini- 
mally edited (for bad traces, resets, etc.). Chi-squareds 
per degree of freedom greater than two were rejected, 
corresponding to a ten percent loss of traces on average. 



VI. MONOPOLE CROSS-SECTIONS AND 
MASS LIMITS 

The analyses of Sec. O above yield a number of "step" 
histograms. Those for set 1, DO aluminum and beryllium, 
are shown in Fig. ^2 for set 2, CDF lead, in Fig.^J and 
for set 3, the CDF aluminum samples, in Fig. ^| after 
taking account of the cuts on the errors discussed below. 

In these plots the variable dS is the standard statisti- 
cal error on the step whereas sS is a systematic error 
which is the rms of the deviations of the actual step from 
the flat fitted steps shown, e.g., in Fig. ^1 If the dis- 
tribution is flat (with little slope) sS is small, whereas 
if it is badly sloped it is large. Bad measurements, e.g., 
very large dipole tails, etc., cannot give a true measure 
of the step and so are removed by this cut. Note that 
Fig. [TUb has a wide dipole response, but still an accept- 
able sS; its wide dipole response comes about because it 
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FIG. 12: Histograms of the steps for the CDF lead samples, 
with an rms dispersion of 0.85 mV. a) 187 layer 4, 5, 6 samples 
along with the 14 remeasurements, having sS < 0.3 mV and 
dS < 0.45 mV, b) 90 layer 2 samples having sS < 0.4 mV 
and dS < 0.6 mV, and c) 95 layer 2 samples having sS < 0.6 
mV and dS < 0.8 mV. Also shown are the four corresponding 
Gaussians for comparison. Note that the distributions in a), 
b), c) are used for limits on magnetic charges of \n\ = 1, 2, 3, 
respectively. 
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FIG. 13: Histogram of the steps from 240 aluminum sam- 
ples (sS < 1.0 mV) from the CDF experiment, with an 
rms dispersion of 2.7 mV (reflecting a higher gain setting of 
SQUID DCIR, see text) compared to the appropriate Gaus- 
sian. These data are used for limits for the three magnetic 
charge values |n| = 2,3, and 6. 



FIG. 14: The error histograms for dS (dashed) and sS (solid) 
for sample sets 1, 2 and 3 (remember the higher gain setting 
for set 3). Note that the the individual (dS,sS) error values 
for the step plots a, b, c, shown in Figure ITUl are (0.22, 0.08), 
(0.34, 0.32) and (0.82, 0.48) mV, respectively. 

is the superposition of two dipoles separated in position. 
Some cases cut by a large sS are due to larger or more 
complicated superpositions of dipole responses. The his- 
tograms of dS and sS are shown in Fig. cuts on dS 
and sS are made in order to control the dispersion of 
the step histograms, and losses of events are taken into 
account in the efficiency e defined below. 

A number of parameters need to be taken from these 
plots and interpreted to yield the cross section and mass 
limits desired. Each plot determines an upper limit 
on the number (N^i) of monopoles of a given magnetic 
charge (n = 1,2,3,6) for each sample set which in turn 
determines an upper limit to the corresponding cross- 
section 



where e is the efficiency for the chosen signal (a monopole 
with charge n) to lie outside a cut excluding smaller val- 
ues of |n|, e is the efficiency of the sample set to cover 
the solid-angle region chosen and to correct for dS and sS 
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Step (mV) 

FIG. 15: A schematic of the upper hmit number analysis pa- 
rameters needed to obtain the cross-section hmits via Tables 
El and EH 



cut losses, and to account for losses of sample material. 
£ is the total luminosity for the p-p exposure delivered 
(172±8pb"^ for DO and 180pb"^ ±5% for CDF 43]). 

The Drell-Yan modelling is that jl^] da/dM^g^ir ^ 
(/3/Mpair)^, where P is the velocity of the monopole in 
the rest frame of the monopole- antimonopole pair. This 
basic cross section, involving the interaction suppression 
factor of is multiplied by the phase space /3 and con- 
volved with the CETQ5 44] quark-parton distributions 
(PDFs). In Ref. 18] we normalized the integral of this 
cross section to (g/e)'^ ~ 4700n^ times the experimental 
CDF and DO pp pijl cross sections ^3 0? ^ind 
obtained the a{MM) Drell-Yan cross section. For Luo's 
thesis an a priori theoretical normalization was made, 
by inserting the quark- ant iquark to mu pairs cross sec- 
tion, 



2J1 



3Q2 



(6.2) 



where is the quark charge and Q is the mass of the 
quark pair, in place of dcj/dM above. The summing over 
quark charges and convolving with PDFs, etc., was also 
done. This procedure agrees with the first one within 
a factor of two, and is stable with respect to choice of 
PDF (e.g., MRST 47]). Equation has integrated 

out a 1 + cos^ center-of-mass angular distribution ap- 
propriate to a jijl virtual photon vertex. The Drell-Yan 
Monte Carlo puts this factor back into the equation and 
throws weighted events according to the mass and an- 
gular distributions of the model. The thrown monopole 
events are Lorentz transformed to the laboratory frame 
and tracked through the various structural materials into 
the appropriate sample (s). We also take into account the 
acceleration or deceleration due to any external magnetic 
fields present (~ 30n GeV/m for CDF). Those events en- 
tering and stopping in the sample layer (s) are summed 
into the "accepted" cross section a a- The acceptance is 
then A = aA/cr, where a is the total cross section without 
any path length or angle cuts. 

In order to extract A/'ui, we must make cuts on the 
step histograms, as shown in Fig. We define a re- 
gion about zero using Co (i.e., ± Co) such that most 



(zero) steps lie there, leaving most of the potential n 7^ 
charges lying outside of this region. The Co to 00 and 
— 00 to —Co region defines truncated Gaussians for mag- 
netic charge ±n, whose standard deviation is assumed 
(in the absence of further information) to be the same as 
that of the central (zero) peak, and shown as the small 
offset Gaussians in the figure. These are centered upon 
the expected positions of the charge n monopoles, which 
is n times the pseudopole calibration value from Table HH 
Sec. mil Then e is the fractional area of these truncated 
Gaussians between Co and infinity. Also, the expected 
number (A^exp) of null events outside ±Co is the frac- 
tion of the central Gaussian lying outside times the total 
number of events in the histogram.^ This number is to 
be compared to the observed number (A/obs) lying out- 
side. We use the Feldman- Cousins tables of A^exp vs TVobs 
to find N^i 48] at 90% confidence level. These parame- 
ters are compiled and displayed in Tables IVl and I VII from 
which (J is calculated and tabulated (the lower limit on 
the monopole mass, mLL, is discussed below). 

The only non-standard cases are the remeasurements 
of sets 1 and 2, in which the loss of possible monopole 
events not remeasured has to be taken into account, by 
another division by £, as remarked in Table IVII We note 
that the cross-section limits range from 0.07 to 9.9 pb, a 
considerable improvement over earlier results from other 
experiments. 

The above results are based on an model using a 
1 + cos^ angular distribution, appropriate for the cou- 
pling of spin-1/2 leptons to photons. As the angular dis- 
tribution of magnetic monopoles is unknown, we try also 
two other distributions to give a fiavor of the changes that 
would occur: isotropic (i.e., 1) and sin^ ^ = 1 — 
These alternative possibilities are given in Table IVIII 
The total Drell-Yan cross-sections for these three angular 
distributions are in the ratio of 1:3/4: 1/2 respectively 
when interpreting as mass limits. 

The cross-section limits above are model dependent, 
but only moderately so, since only the shape of the 
dcr/dMpair distributions is relevant to the acceptance 
A, apart from the straightforward angular distribution, 
dE/dx, and magnetic field tracking considerations. The 
interpretation of these limits as monopole mass limits 
is directly model dependent. We have used the Drell- 
Yan model for lack of any field theoretic results. (For 
the status of the theory, see Ref. 0.) We modify the 
couplings from /Tic to jTic (as appropriate in dE / dx 
calculations) rather naively, but we have also included 13^ 
velocity suppression factors, a conservative choice. More- 
over, we note that there is a unitarity limit that comes 
in at n ~ 3, and we use n = 3 cross-sections as the uni- 
tarity limit for all n > 3 (i.e., we use three cross-sections, 
those for n = 1 and n = 2, and that for n = 3 for higher 



^ The value of Co is chosen to minimize A^exp while retaining good 
efficiency for finding n 7^ 0. 
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TABLE V: Summary of parameters for the various sample sets and the various magnetic charge values (n) of the step histograms 
needed via Fig. 1151 Here "response" refers to the mean value expected for a monopole of strength n. The table lists only those 
samples that have sufficient sensitivity for that value of n. 



Set 


n 


Dispersion 


Co 


Response 




A^exp 


Nobs 






(mV) 


(mV) 


(mV) 








1 Al 


1 


0.73 


1.47 


2.46 


222 


9.8 


8 


1 Al RM 


1 


0.73 


1.47 


2.46 


8 


0.4 





2 Pb 


1 


0.85 


1.6 


2.46 


187 


11.2 


14 


2 Pb RM 


1 


0.85 


1.6 


2.36 


14 


0.8 





1 Al 


2 


0.73 


3.0 


4.92 


222 








2 Pb 


2 


0.85 


3.5 


4.92 


90 








3 Al 


2 


2.7 


8.8 


10.64 


240 


0.3 





1 Be 


3 


0.73 


3.0 


7.38 


6 








2 Pb 


3 


0.85 


3.5 


7.38 


95 








3 Al 


3 


2.7 


8.8 


16.0 


240 


0.3 





1 Be 


6 


0.73 


3.0 


14.8 


6 








3 Al 


6 


2.7 


8.8 


31.9 


240 


0.3 






TABLE VL The remaining parameters and final limits (90% confidence level) from Drell-Yan modeling for cross-sections and 
monopole masses. The upper limit A^ui comes from a Feldman- Cousins analysis 48]. Note: Set I Al RM and 2 Pb RM (lines 2 
and 4) are corrected for losses of samples not remeasured by dividing by e again. 



Set 


n 


iVul 

(90% CL) 


e 


e 


^+1 


C 

(Pb-^) 


-^1 
(pb) 


m\\ 
(GeV/c") 


1 Al 


1 


4.4 


0.91 


0.94 


0.026 


172 


1.2 


250 


1 Al RM 


1 


2.0 


0.91 


0.94 


0.026 


172 


0.6 


275 


2 Pb 


1 


10.3 


0.84 


0.53 


0.013 


180 


9.9 


180 


2 Pb RM 


1 


1.7 


0.82 


0.53 


0.011 


180 


2.4 


225 


1 Al 


2 


2.4 


1.00 


0.94 


0.007 


172 


2.1 


280 


2 Pb 


2 


2.4 


0.95 


0.83 


0.017 


180 


1.0 


305 


3 Al 


2 


2.1 


0.75 


0.88 


0.10 


180 


0.2 


365 


1 Be 


3 


2.4 


1.00 


1.0 


0.0036 


172 


3.9 


285 


2 Pb 


3 


2.4 


1.00 


0.87 


0.028 


180 


0.5 


350 


3 Al 


3 


2.1 


0.997 


0.88 


0.19 


180 


0.07 


420 


1 Be 


6 


2.4 


1.00 


1.0 


0.013 


172 


1.1 


330 


3 Al 


6 


2.1 


1.00 


0.88 


0.057 


180 


0.2 


380 



TABLE VIL Alternative interpretations for different production angular distributions of the monopoles, comparing 1 and 
1 — cos^ to the (repeated here) 1 + cos^ limits. Here the acceptance Aa corresponds to the distribution 1 + acos^ 0, and 
similarly for the cross section and mass limits (all at 90% confidence level). 



Set 


n 




LL 

m+i 


Ao 




LL 

mo 


A-i 










(pb) 


(GeV/c") 




(pb) 


(GeV/c^) 




(pb) 


(GeV/c") 


1 Al 


1 


1.2 


250 


0.024 


1.2 


240 


0.021 


1.4 


220 


1 Al RM 


1 


0.6 


275 


0.024 


0.6 


265 


0.021 


0.7 


245 


2 Pb 


1 


9.9 


180 


0.011 


12 


165 


0.0055 


23 


135 


2 Pb RM 


1 


2.4 


225 


0.009 


2.9 


210 


0.0045 


5.9 


175 


1 Al 


2 


2.1 


280 


0.0068 


2.2 


270 


0.0060 


2.5 


250 


2 Pb 


2 


1.0 


305 


0.018 


0.9 


295 


0.016 


1.1 


280 


3 Al 


2 


0.2 


365 


0.10 


0.2 


355 


0.096 


0.2 


340 


1 Be 


3 


3.9 


285 


0.0025 


5.6 


265 


0.0003 


47 


180 


2 Pb 


3 


0.5 


350 


0.029 


0.5 


345 


0.031 


0.5 


330 


3 Al 


3 


0.07 


420 


0.20 


0.07 


410 


0.24 


0.06 


405 


1 Be 


6 


1.1 


330 


0.008 


1.7 


305 


0.0008 


18 


210 


3 Al 


6 


0.2 


380 


0.066 


0.2 


375 


0.082 


0.2 


370 
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FIG. 16: The curves are Drell-Yan cross sections versus 
monopole mass with cross-section upper hmits (90% CL) in- 
terpreted as mass hmits (cross section upper hmit intersects 
Dreh-Yan curve at estimated lower limit of mass as shown 
by + and • markers). Three possible center-of-mass angular 
distributions are considered, of form 1 + acos^ 0, with a = 1, 
and —1, respectively. 



n). Then, taking the intersection of the a limits with 
the Drell-Yan curves we obtain the lower limit for the 
mass (to the nearest 5 GeV) of each given monopole case, 
as entered in Tables IVll and IVlll and shown in Fig. ^1 
These mass limits are an improvement of a factor of about 
two or more over previous direct experimental results. 



VII. SUMMARY 

We have improved the cross-section and mass limits on 
directly produced magnetic monopoles over the best pre- 
vious other limits of Bertani et al. 13] and the lunar rock 
limits of Ross et al. 17]. We have improved and extended 
our own prior limits from DO samples 18] to those from 
CDF, with a consistent overall analysis, presented in this 
paper. 



Since limits are given, it is hard to assign errors. How- 
ever, we can give a range or "spread" of the values from 
the quoted means to give a flavor of the variations to be 
expected. These spreads are in general a combination of 
statistical and systematic errors. Since we are modelling 
according to the Drell-Yan process, the systematic errors 
on the acceptances are unknown; it would take the eval- 
uation of another reasonable proposed model, or a true 
quantum field theoretic prediction, to assess that. For ex- 
ample, a factor of two change in the cross section would 
lead to a change in the mass limits of about 20 GeV. The 
acceptances are estimates in any case, since the modeling 
of the detector elements, the dE/dx losses, and the track- 
ing of the particles through the various elements are not 
exact. The statistical error on the acceptances are 3% in- 
dividually and 1% from smoothing over a mass interval 
(with a quadratic functional fit) of some 9-11 calculations 
spanning a 200-300 GeV mass region. The efficiencies {e 
and e) are good to about 2% and the luminosities to some 
5%. Combining 3%, 2%, 2%, and 5% quadratically gives 
a spread of ~ ±7% on in addition to the variations 
due to differing assumed angular distributions. Since the 
Dress- Yan cross sections are closely linear on a semilog 
plot against mass, this converts to a ~ ±1 GeV spread in 
mass. The unknown angular distributions give a greater 
uncertainty. 

Since the true angular distributions are unknown, we 
take as a "base" an isotropic one (at 90% CL). In addi- 
tion, we choose the largest mass limit for each value of 
the magnetic charge n. These come from DO Al set 1 
for n = 1, and CDF Al set 3 for n = 2, 3 and 6. From 
Fig.mand TablelVTll we see that these have cross section 
limits of 0.6, 0.2, 0.07, and 0.2 pb, respectively. The cor- 
responding mass limits are 265, 355, 410, and 375 GeV, 
with spreads of (+10,-20), (+10,-15), (+10,-5), and 
±5, respectively, due to replacing isotropic distributions 
by ones of the form 1 ± cos^ 6. 

These cross section limits are some 250 to 2500 times 
smaller and the mass limits are about 3 times larger than 
those of Bertani et al. |13l|. a nd a significant improvement 
over those of Ross et al. jlTl] . as well. 
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APPENDIX: SIMPLIFIED THEORY OF 
MONOPOLE DETECTOR 

This appendix describes the basis of the functioning 
of our magnetic monopole dectector. It works by detect- 
ing the magnetic flux intercepted by a superconducting 
loop contained within a superconducting cylinder. The 
detector is sketched in Fig. El 

In order to incorporate finite-size effects, we consider 
first a perfectly conducting right circular cylinder of ra- 
dius a of semi-infinite length, with axis along the z-axis, 
and with a perfectly conducting circular bottom cap at 
z = 0. We use cylindrical coordinates p, 0, and z. 

Because the boundaries are superconductors, the nor- 
mal component of B must vanish on the surfaces, that 
is, 



Br 



0, B, 



(A.l) 



z=0 



Now suppose a magnetic pole of strength g is placed 
on the z axis dit z = z' > 0. This could either be a mag- 
netic monopole (magnetic charge) or one pole of a very 
long electromagnet ("pseudopole"). Imagine a circular 
conducting loop of radius r < a centered on the axis of 
the cylinder and perpendicular to that axis, with center 
dit z = Z. Inside the cylinder and outside of the loop B 



is derivable from a magnetic scalar potential, 

B = -V(/)M, (A.2) 

since we may ignore the displacement current, because 
the time variation is negligible. 0m satisfies Poisson's 
equation, in cylindrical coordinates: 



V B 



1 d d 1 



4:7rgd{T — r'). 



(A.3) 



where is the position of the monopole, = (p', 0\ z'). 
This is the equation for a Green's function, which we can 
express in separated variables form. That is, we write 



DM = — I dk cos kz cos kz' 

Jo 



E ^e^"'^'-''^9m{p,p';k), (A.4) 



m= — oo 



where, in view of the first boundary condition in 
Eq. (lA.lj) . we may express the reduced Green's function 
in terms of modified Bessel functions: 

gm{p,p']k) = -47rglm{kp<) 

Km{kp>) - Im{kp>)^§^ 

(A.5) 

where p< (p>) is the lesser (greater) of p, p' . If the 
monopole is confined to the z axis, only the m = term 
survives: 



6m = / dk cos kz cos kz' 

Jo 



Ko{kp) ^ Io{kp) 



Ki{ka) 
h{ka) 



(A.6) 



which uses 

I'^{x)=h{x), K'^{x) = -Ki{x). (A.7) 
By integrating over the cross section of the loop using 



Jo 



dttKo{t) = -xKi{x) + 1, (A.8a) 
dttIo{t) = xli{x), (A.8b) 

we obtain the following formula for the magnetic flux 
subtended by the loop, 

^ = J dS-B = 47Tg [r]{Z - z') - F{Z, z')] , (A.9) 

where the step function is 
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and the response function is 

^/ /X 2r ^ , z z' 

Fiz^z) = / dxsmx—cosx — 

7T a Jq a a 

X ^^Kl{xr/a) - ' 

(A.ll) 

Now suppose that the pole is slowly moved from a 
point far above the loop, z' = +00, to a point below the 
loop, z' = Zq^ Z > zq. Then from Maxwell's equation 

Id 47r 
VxE = — -B-— J^, (A.12) 

cot c 

where Jm is the magnetic current density, the emf in- 
duced in the loop is 

/d^ 47r 
^■dl,= -— + —g5{t), (A.13) 

if t = is the time at which the pole passes through the 
plane of the loop. The net change in emf gives rise to a 
persistent current / in the 5iipe?iconducting loop, 

1 47r 47r 

LI= / £dt = —A^^—g = —gF{Z,zolA.U) 
J-00 c c c 

where L is the inductance of the loop, and the response 
function F is given in Eq. (|A.11|) . This is just a statement 
of the Meissner effect, that the flux change caused by the 
moving monopole is cancelled by that due to the current 
set up in the loop. 

Now suppose that the loop is very far from the bottom 
cap, Z a. Then only small x contributes to the integral 
in Eq. (|A.11|) . and it is easy to see that 



consider, instead, an infinite, open-ended, superconduct- 
ing cylinder, with the monopole passing from z = +00 to 
z = —00, at either extreme half the flux must cross the 
plane of the loop, so with the uniformity assumption we 
get the same result: 



LI = 



47rg 



c 



Aug 



c \ 



(A.17) 



The simple assumption of a uniform magnetic field is 
apparently justified by the exact result (|A.15|) . 

We conclude this appendix by noting how the exact 
calculation is modified for an infinite superconducting 
cylinder. In the magnetic scalar potential, the integral 
over k mode functions in Eq. (|A.4|) is replaced by 



ik{z — z') 



(A.I8) 



which has the effect of replacing the flux expression ()A.9|) 

by 



<^> = 2ng[e{z-z')-F{Z-z',0)], 



where 



e{x) 



r 1, X 

1-1, X 



< 0. 



(A.19) 



(A.20) 



Then the induced current in the detection loop when the 
monopole passes from a point above the loop z' = Z -\- ^ 
to a point, equidistant, below the loop, = Z — ^, is 



Ang 



Ang 



LI=^F{^,0)^^(1-^ 



(A.21) 



£dt 



Airg 



a2 



(A.15) 



SO the signal is maximized by making the loop as small 
as possible, relative to the radius of the cylinder. We get 
the full flux of the monopole only for a loop in empty 
space, a/r 00. This perhaps counterintuitive effect is 
due to the fact that the superconducting walls confine 
the magnetic flux to the interior of the cylinder. Thus 
for the superconducting can, the induced current in the 
detection loop caused by the passage of a monopole from 



00 to = is 



LI = 



Aug A<l> Aug ^{z' = 0) 



(A.16) 



where the last limit applies if ^/a ^ 1. This result co- 
incides with that in Eq. (|A.15|) . The function R{^) = 
|F(^,0)/(1 — r'^/a'^) + |, corresponding to a monopole 
starting from a point zi far above the loop, — Z a, 
and ending at a point = Z — ^, is plotted as a func- 
tion of ^ for our parameter values in Fig. [7t, where it 
is shown to agree well with experimental data. This re- 
sponse function coincides with the result obtained from 
Eq. ()A.14|) . because 



-F(2Z 
2 ^ 



2\ 1 



F(e,0), (A.22) 



which yields the result (|A.15|) if one assumes that the 
magnetic field is uniform across the can's cross section at 
the position of the loop when the pole is at the bottom, 
because all the flux must pass up through the can. If we 



if Z/a ^ 1. This shows that the effect of the endcap 
(which of course is not present in actual detector) is neg- 
ligible, demonstrating that the fact that the supercon- 
ducting shield is of finite length is of no significance. 
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